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Abstract 

To understand the role of electron correlation and disorder on persistent currents in isolated ID rings 
threaded by magnetic flux </>, we study the behavior of persistent currents in aperiodic and ordered binary 
alloy rings. These systems may be regarded as disordered systems with well-defined long-range order so 
that we do not have to perform any configuration averaging of the physical quantities. We see that in the 
absence of interaction, disorder suppresses persistent currents by orders of magnitude and also removes its 
discontinuity as a function of </>. As we introduce electron correlation, we get enhancement of the currents 
in certain disordered rings. Quite interestingly we observe that in some cases, electron correlation produces 
kink-like structures in the persistent current as a function of (jj. This may be considered as anomalous 
Aharonov-Bohm oscillations of the persistent current and recent experimental observations support such 
oscillations. We find that the persistent current converges with the size of the rings. 
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1 Introduction 

In the mesoscopic and nanoscopic world phase co- 
herence of the electronic states are of fundamental 
importance, and the phenomenon of persistent cur- 
rent is a spectacular consequence of quantum phase 
coherence in these regimes. In a pioneering work 
Biittiker, Imry, and Landauer p suggested that 
even in the presence of impurity, a small conduct- 
ing isolated ring enclosing a magnetic flux cj) car- 
ries a current in the ground state, a current which 
persists (does not decay) in time, and periodic in 
(j) with periodicity </>o — ch/e^ the elementary flux 
quantum. Since then the phenomenon of persis- 
tent current in mesoscopic systems has been dis- 
cussed quite extensively in the literature both the- 
oretically EiEiiiisiEiiziiHiEiiiniiniiiiiisiiiiini 

as well as experimentally [TBI El UHl [HOni ■ How- 
ever, till now the observed experimental features of 
the persistent currents are not well-understood the- 
oretically. A typical example of such discrepancy 
between theory and experiment is that the ampli- 
tude of the measured persistent currents are orders 
of magnitude larger than the theoretical predictions. 
It is believed that the electron-electron correlation 
and disorder have major role on the enhancement of 
persistent currents, but no consensus has yet been 
reached. Another controversial issue is that experi- 
mentally both 00 a-iid 0o/2 periodicity has been ob- 
served, and it is also found that the '/)o/2 oscilla- 
tions near zero magnetic field exhibit diamagnetic 
response. The explanation of these results in terms 
of the ensemble averaged persistent currents is also 
quite intriguing, and the calculations show that the 
disorder averaged current crucially depends on the 
choice of the ensemble 01 El • In order to reveal the 
role of disorder and electron correlation on the per- 
sistent currents, in this work we focus attention on 
certain systems which closely resemble the disorder 
systems where we do not require any configuration 
averaging. These are chemically modulated struc- 
tures possessing well-defined long-range order, and, 
as specific examples we consider the aperiodic and 
ordered binary alloy rings. We confine ourselves to 
small ID rings where persistent currents can be cal- 
culated exactly, and we obtain many interesting new 



results as a consequence of electron correlation and 
disorder. One such result is the enhancement of per- 
sistent currents in these systems due to electron cor- 
relation. Another important observation is the ap- 
pearance of kink-like structures in the persistent cur- 
rent due to electron-electron correlation, which may 
be considered as anomalous Aharonov-Bohm oscilla- 
tions. With the recent advancements in sub-micron 
technology, such systems can be easily fabricated 
in the laboratory, and in fact, in a recent experi- 
ment Keyser et al. |2U| reported similar anomalous 
Aharonov-Bohm oscillations from the transport mea- 
surements on small rings with less than ten electrons. 
Our study may also be helpful to understand the 
physical properties of Benzene-like rings, and other 
aromatic compounds in the presence of magnetic flux. 

This paper is organized as follows. In section 2, we 
present the calculation of persistent currents in or- 
dered binary alloy rings and investigate their behav- 
ior in the presence of Coulomb repulsion. In section 
3, we describe our results for the incommensurate 
rings in the presence of electron-electron interaction. 
Lastly, we conclude in section 4. 

2 Ordered binary alloy rings 

In this section we describe the current-flux charac- 
teristics for the ordered binary alloy rings at T = 0. 
We use the tight-binding Hubbard Hamiltonian with 
a pure Aharonov-Bohm flux (p (in units of 4>q) with- 
out any Zeeman term. The magnetic vector potential 
modifies the hopping integral by a phase factor and 
the Hamiltonian for a TV-site ring becomes 

N N 

H = <^^(4^c,+i,„e'®-h/i.c.) + C/^?i,|n,; 

N-1 

+ (eA"j,a- + eB7ij+l,a) (1) 

cr 1=1,3,... 

Here c\^{cia) is the creation (annihilation) operator 
and riicy is the number operator for the electron in the 
Wanniar state \ia >. The parameters t and U are re- 
spectively the nearest-neighbor hopping integral and 
the strength of Hubbard correlation, and, ca and 
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Figure 1: Schematic diagram of a ID ordered binary 
alloy ring threaded by magnetic flux (j). 

are the site potentials for the A and B type atoms. 
The phase factor is = 2'K(j)/N . Henceforth, we take 
t — —1 and use the units c ^ e = h = I. We always 
choose N to be even so as to preserve the perfect bi- 
nary ordering of the two types of atoms in the ring 
(see Fig-P). 

At zero temperature, persistent current in an iso- 
lated ring threaded by magnetic flux is given by 



(2) 



where Eo{(p) is the ground state energy. We calcu- 
late I{(f>) exactly by numerical diagonalization of the 
Hamiltonian. 

Let us now study the behavior of persistent cur- 
rents in the ordered binary alloy rings, and investi- 
gate the role of electron-electron interaction on the 
currents. In a pure ring consisting of either ^ or _B 
type atoms without electron correlation, the persis- 
tent current as a function of 4> is discontinuous at cer- 
tain points due to ground state degeneracy, and, the 
I — (/) curve exhibits a saw-tooth like behavior[3]. In a 
binary alloy ring with U = 0, this discontinuity com- 
pletely disappears as illustrated in Fig. Ela) by the 
solid curve. The reason for this is that the binary al- 
loy configuration may be considered as a perturbation 
over the perfect ring which lifts the ground state de- 
generacy, and consequently, makes /((/>) a continuous 





Figure 2: The / — curves of four (T, T, i) electron 
ordered binary alloy rings, (a) = 4, TVg = 4. The 
solid, dotted, small dashed and dashed lines corre- 
spond to U — 0, 2, 4 and 10 respectively, and (b) 
N — 8, Ae = 4. The solid, dotted, small dashed and 
dashed lines are respectively for U = 0, 2, 4 and 6. 



function of (p. As we switch on the electron-electron 
interaction in a pure ring, persistent current always 
decreases with the increase of U. However, in the 
ordered binary alloy rings, depending on the number 
of electrons A^e, we observe enhancement of the per- 
sistent current for low values of but eventually it 
decreases when U becomes very large. Such a situ- 
ation is depicted in Fig. Efa) where we display the 
I — (j) curve for a half-filled ordered binary alloy ring 
with four electrons (two up and two down spin elec- 
trons). Here we see that for U = 2 (dotted line) and 
U — 4 (small dashed line), the current amplitudes 
are significantly larger than the non-interacting case, 
whereas for U =10 (dashed line) the current am- 
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plitude is less than that from the U = case. This 
enhancement takes place above quarter-filling, i.e., 
when Ne > N/2, and can be easily understood as 
follows. As N is even, there are exactly N/2 number 
of sites with the lower site potential energy. If we 
do not take into account the electron-electron inter- 
action then above quarter-filling, it is preferred that 
some of these lower energy sites will be doubly occu- 
pied in the ground state. As we switch on the Hub- 
bard correlation, the two electrons that are on the 
same site repel each other and thus causes enhance- 
ment of persistent current. But for large enough U, 
the electron-hopping is strongly suppressed by inter- 
action, and the current gets reduced. On the other 
hand, at quarter-filling and also below quarter-filling, 
no lower energy site will be doubly occupied in the 
ground state and hence there is no possibility of en- 
hancement of persistent current due to Coulomb re- 
pulsion. In these systems we always get suppression 
of persistent current with the increase of interaction 
strength U. In Fig.|2Ib), we plot the I—<j) curves for a 
quarter-filled binary-alloy ring with U = 0, 2, 4 and 
6, and it clearly shows that the persistent currents 
are always suppressed by interaction. 

Let us now describe the behavior of persistent cur- 
rent with system size N in the ordered binary-alloy 
rings keeping N^/N constant. For this purpose we 
calculate current amplitude /q at some typical value 
of magnetic flux = 0.25 and in Fig.|3|we plot the /q 
versus N curves. The results for the non-interacting 
rings are presented in Fig. |3Ja) , where the solid and 
dashed lines correspond the rings of size N — AN^ 
and — ANe + 2 respectively. Since the dimension 
of the Hamiltonian matrices for the interacting sys- 
tems increases very sharply with A^ for higher num- 
ber of electrons A'e and also the computational op- 
erations are so time consuming, we present the vari- 
ations for the interacting rings with size A^ = 2Ne 
only. In Fig.lSJb) the results for the interacting rings 
are plotted. The solid and dotted lines, representing 
the results respectively for the rings with U = 0.2 
and U = 0.6, almost coincide with each other, while, 
those results for the rings with U — 1.5 and U — 3 
are respectively represented by the small dashed and 
dashed lines. These results can predict the variations 
of current amplitude for the rings with size A^ = ANe 




Figure 3: Current amplitude Iq as a function of sys- 
tem size N in ordered binary- alloy rings, (a) The 
solid and dashed lines correspond the non-interacting 
rings of size A^ = 4A"e and A^ = 4A'e -I- 2 respectively, 
(b) The solid, dotted, small dashed and dashed lines 
correspond the interacting rings of size A^ = 2Ae with 
U = 0.2, 0.6, 1.5 and 3 respectively. 

and also A^ — 4A'e + 2. It is apparent from Fig. O that 
the current amplitude gradually decreases with sys- 
tem size i.e., we get a converging behavior of current 
amplitude with A^ and most interestingly we see that 
in the interacting rings current amplitude converges 
to zero for any non-zero value of U. These results 
predict that in a realistic bulk system Iq goes to zero 
as soon as the interaction is turned on. 

3 Rings with incommensurate 
site potentials 

In this section, we study persistent currents in ID 
rings with quasi-periodic site potentials, and investi- 
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gate the effects of electron-electron interaction on the 
currents. We describe a iV-site ring with incommen- 
surate site potentials by the Hamiltonian 

N 

<T i=l 

N N 

+ ^ ^ e cos(iA7r)4^^Ci,^ + ^ X! "»T"-ii (3) 

where A is an irrational number, and as a typical ex- 
ample we take it as the golden mean ^ ^"^2^ ) • Setting 
A = we get back the pure ring with identical site 
potential e. 

To understand the precise role of electron-electron 
interaction on persistent current in the presence of 
incommensurate site potentials, let us first consider 
the simplest possible system which is the case of a 
ring with two opposite spin electrons. In Fig. ^ we 
plot the I — (j) curves for a 30-site incommensurate 
ring with U = (solid line), U — 1 (dotted line), and 
[/ = 3 (dashed line). In the absence of any inter- 
action, persistent current is greatly reduced by the 
incommensurate site potential, and in fact. Fig. ^ 
shows that the I — 4> curve for the non-interacting 
case (solid line) almost coincides with the abscissa. 
This result can be easily understood from the argu- 
ment that in the presence of aperiodic site poten- 
tials, the electronic eigenstates are critical |^ 122 
which tends to localize the electrons, and thus re- 
duces the current. But this situation changes quite 
dramatically as we switch on the electron-electron in- 
teraction. Fig. 01 clearly shows that electron correla- 
tion considerably enhances persistent current for low 
values of U. This is because the repulsive Coulomb 
interaction does not favor double occupancy of the 
sites in the ground state, and also it opposes con- 
finement of the electrons due to localization. Thus 
the mobility of the electrons increases as we intro- 
duce Hubbard correlation and gives enhancement of 
persistent current. But this enhancement ceases to 
occur after certain value of U due to the ring geom- 
etry, and persistent current then decreases as we in- 
crease U further. We also observe that some strange 
kink- like structures appear in the / — </> characteristics 
around (j) — ±0.5 for non-zero values of U. Quite sur- 
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Figure 4: The I — 4> curves for two (T, i) electron in- 
commensurate rings with N — 30. The solid, dotted 
and dashed lines are respectively for [/ = 0, 1, and 3. 



prisingly we notice that the persistent currents inside 
these kinks are independent of the strength of Hub- 
bard correlation U . Let us now analyze this result. 
For two opposite spin electrons, total spin S can have 
the values S* = and 1 . The Hamiltonian of this sys- 
tem for any (j) can be block diagonalized by proper 
choice of the basis, and this can be achieved by tak- 
ing all the basis states in one sub-space with 5 = 0, 
while those in the other sub-space with S ~ 1. It 
is easy to see that in the sub-space spanned by the 
basis set with 5* = 1, the block Hamiltonian is free 
from J7, and, the corresponding energy eigenvalues 
and eigenstates are [/-independent. In the absence 
of interaction, these [/-independent energy levels are 
always above the ground level for any 0. But for non- 
zero values of [/, one of these [/-independent energy 
levels becomes the ground state energy of the system 
in certain domains of (p. In these regions we have 
kinks in the I — (p curves, and it is obvious that the 
persistent currents inside these kinks are independent 
of the Hubbard correlation U. We observe that inter- 
action does not alter the (f>o periodicity of persistent 
current in this two-electron system. 

Next we consider incommensurate rings with two 
up and one down spin electrons as the representa- 
tive examples of three-electron systems. The I — (j) 
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characteristics for the half-filled {i.e., N = 3, ~ 3) 
system with [7 = 4 are shown in Fig.|2fa). In this fig- 
ure, the dotted line corresponds to a pure ring which 
exhibits discontinuous jumps at = 0, ±0.5 due to 
the crossing of the energy levels. Interestingly, we ob- 
serve that this pure (A = 0) half-filled three-electron 
system exhibits a perfect (/'o/2 periodicity, and we 
will see that this is a characteristic feature of the 
pure half-filled systems with odd number of electrons. 
It is evident from the solid curve of Fig. EJa) that 
this 00 /2 periodicity no longer exists as we introduce 
the incommensurate site potentials, and we have the 
usual (f>o periodicity. Moreover, in this case I{cf>) be- 
comes a continuous function of (j) as the perturbation 
due to disorder lifts the ground state degeneracy at 
the crossing points of the energy levels. The charac- 
teristic features of the persistent currents are quite 
different in the non-half-filled rings with two up and 
one down spin electrons, and in Fig.[5^b) we present 
the results for a 12-site ring with incommensurate 
site potentials. The solid line is the I — <j) curve for 
the non-interacting electrons {U — 0), while the dot- 
ted and dashed lines represent the I — 4> curves for 
U — 4: and U = 50 respectively. The role of electron- 
electron interaction on persistent current in the pres- 
ence of incommensurate site potentials becomes evi- 
dent from these curves. For low values of U, the / — (/) 
curve resembles to that for the non-interacting case, 
and the persistent currents do not show any discon- 
tinuity. But for large enough U, kink-like structures 
appear in the I — cf) characteristics as illustrated in 
Fig. Elb) by the dashed line, e.g., around the point 
= 0. In this case also, the kinks are due to the 
[/-independent eigenstates like the two electron sys- 
tems, and as explained earlier, the currents inside the 
kinks are independent of the strength of correlation 
U . It is observed that the persistent currents always 
have the 0o periodicity in the non-half-filled systems. 
We also notice that for the half-filled systems, persis- 
tent currents always decrease as we increase C/, while 
in the non-half-filled rings we obtain significant en- 
hancement of the currents due to interplay between 
electron correlation and the incommensurate site po- 
tentials. 

As a systematic approach, next we investigate the 
behavior of persistent currents in four-electron sys- 
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Figure 5: The I — (p curves for three (t,Tji) elec- 
tron incommensurate rings, (a) Half-filled {N = 3) 
systems with U = 4. The dotted and solid lines cor- 
responds to A = and A ^ respectively, (b) Non- 
half-filled {N = 12) systems with A ^ 0. The sohd, 
dotted and dashed lines are respectively for [/ = 0, 4, 
and 50. 



tems with incommensurate site potentials, and as the 
representative examples we consider rings with two 
up and two down spin electrons. The I — (p curves 
for the half-filled systems are plotted in the Fig.lHta). 
The solid, dotted and dashed lines are for the cases 
with U = 0, U = 4 and U — 10 respectively. It 
is evident from these curves that the current ampli- 
tude gradually decreases with the increase of interac- 
tion strength U. This indicates that in the half-filled 
systems, the mobility of the electrons gradually de- 
creases with the increase of U, and we see that for 
large enough U, the half-filled system goes to an in- 
sulating phase. This kind of behavior holds true in 
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Figure 6: The I — cf) curves for four (T,Tiiii) in- 
commensurate rings, (a) Half-filled {N = 4) sys- 
tems. The solid, dotted and dashed lines corresponds 
[/ = 0,4 and 10 respectively, (b) Non-half-filled 
{N = 8) systems. The solid and dashed lines are 
respectively for U — A and 16. 



any half- filled system because at large enough U, ev- 
ery site will be occupied by a single electron and the 
hopping of the electrons will not be favored due to 
strong electron-electron repulsion. In Fig. E^b), we 
display the I — (j) curves for the non-half-filled four- 
electron systems with the aperiodic Harper potential. 
The solid and dotted lines are the I — (p curves for a 
8-site ring with U = A and C/ = 10 respectively. This 
figure depicts that for low U persistent current 
has no discontinuity, but kinks appear in the I — (f> 
curve at large U. These kinks arise at sufficiently 
large U due to additional crossing of the ground state 
energy levels as we vary (f). It may be noted that in 
the present case kinks are due to the {/-dependent 



states, and not from the iJ-independent states as in 
the previous cases. Both for the half-filled or non- 
half-filled incommensurate rings with four electrons, 
we observe that the persistent currents always exhibit 
4>Q periodicity. 
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Figure 7: The / — curves for five (T, T; T: i) elec- 
tron incommensurate rings. The curves are for the 
non-half-filled {N = 7) systems. The solid and dot- 
ted lines are respectively for U = 18 and 120. 

Lastly, we consider five electron aperiodic rings, 
and calculate persistent currents in rings with three 
up and two down spin electrons. In the pure half- 
filled case {i.e., N ^ 5, = 5 and A = 0) we get 
4>o/'2 periodicity in persistent current, and we have 
already observed such period halving in other pure 
half-filled systems with odd number of electrons {e.g., 
N = 3, Ne = 3 and A = 0). Like the three-electron 
half-filled incommensurate systems, also in this case 
the (f>o periodicity of persistent current is restored 
once we introduced incommensurate site potentials. 
The I — (j) curves for the non-half-filled five-electron 
systems with N — 7 are shown in Fig. [T] The solid 
and dotted lines are respectively for the cases with 

= 18 and J7 = 120. As in the non-half-filled three- 
electron system, here also kinks appear in the I — 4> 
curves above some critical value of U . Another im- 
portant observation is that for large U {U = 120), 
the maximum amplitude of the current remains finite. 
This is quite natural since we are considering the sys- 
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tems with N > where some sites are always empty 
so that electrons can hop to the empty site, and thus 
the system always remains in the conducting phase. 
We also see that in these non-half-filled five-electron 
rings persistent currents always have the 4>q period- 
icity. 
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Figure 8: Iq versus N curves for the systems with 
incommensurate site potentials, (a) Non-interacting 
rings, where the solid and dashed lines correspond 
the rings of size N = 2Ne and N = ANe respectively, 
(b) The solid, dotted, small dashed and dashed lines 
correspond the interacting rings of size N = 2Ne with 
U = 0.2, 0.6, 1.5 and 6 respectively. 

Now we address the behavior of persistent current 
with system size N keeping N^/N constant, for the 
rings with incommensurate site potentials. Like or- 
dered binary-alloy rings here we determine the cur- 
rent amplitude /q at some typical value oi (j) — 0.25 
and Fig. |S1 displays the Iq versus N curves. The 
solid and dashed lines in Fig.lSJa) respectively corre- 



spond the results for the non-interacting rings of size 
N = 2Ne and N = AN^. On the other hand the vari- 
ations of the interacting rings with size N = 2Ne are 
shown in Fig. IHI^b). In this figure the sofid, dotted, 
small dashed and dashed lines correspond the results 
for the rings with U = 0.2, 0.6, 1.5 and U = 6 re- 
spectively. These curves reveal that the current am- 
plitude gradually decreases with system size N and 
for large value of TV it eventually drops to zero value. 
Thus here we also get a converging nature of current 
with the size of the systems and we can say that the 
current amplitude /o converges to zero for any non- 
zero value of U for large N. Thus these results also 
predict that in a realistic bulk system Iq vanishes as 
soon as the interaction is switched on. 

4 Conclusions 

In conclusion, we have studied exactly the charac- 
teristic features of persistent currents in aperiodic 
and ordered binary alloy rings in the presence of 
electron-electron interaction. These systems, which 
are neither pure nor disordered rings, have well- 
defined structural order, and, our study yields many 
interesting results due to interplay between electron- 
electron interaction and disorder in these systems. 
Our main results are : i) In the absence of inter- 
action, the discontinuity in /(</>) as a function of 
disappears due to disorder. This has been observed 
both in the ordered binary alloy rings and also in the 
aperiodic rings, ii) In pure rings with electron cor- 
relation, we observe both <f>Q and 4'o/2 periodicities 
in the persistent currents. However, in the incom- 
mensurate and ordered binary alloy rings persistent 
currents always have the periodicity. Hi) In the or- 
dered binary alloy rings, above quarter-filling we get 
enhancement of persistent current for small values of 
[/, and the current eventually decreases when U be- 
comes large. On the other hand, below quarter-filling 
and also at quarter-filling, persistent current always 
decreases with the increase of U . iv) An important 
finding is the appearance of kink- like structures in the 
I — (j) curves of the incommensurate rings only when 
we take into account electron-electron interaction. 
Quite surprisingly we observe that in some cases the 



8 



persistent currents inside the kinks are independent 
of the strength of the interaction. These kinks give 
rise to anomalous Aharonov-Bohm oscillations in the 
persistent current, and recently Keyser et al. [2Dj ex- 
perimentally observed similar anomalous Aharonov- 
Bohm oscillations in the transport measurements on 
small rings, v) The current amplitude gradually de- 
creases with N both for the non-interacting and in- 
teracting rings, keeping Ne/N constant, i.e., we get 
a converging behavior of current with the size of the 
rings. Most interestingly, we predict that in the re- 
alistic bulk systems Iq vanishes as soon as the inter- 
action is switched on. Thus the current amplitude 
vanishes for a ring of macroscopic size. 
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